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Abstract. A 50 min time series of one-dimensional slit- spectrograms, taken in quiet sun at disk centre, observed at the German 
, Vacuum Tower Telescope (Observatorio del Teide), was used to study the global and spatial variations of different line pa- 

■ rameters. In order to determine the vertical structure of the photosphere two lines with well separated formation heights have 
been considered. The data have been filtered of p-modes to isolate the pure convective phenomenon. From our studies of global 
correlation coefficients and coherence and phase shift analyses between the several line parameters, the following results can 
be reported. The convective velocity pattern preserves structures larger than 1 "0 up to the highest layers of the photosphere (~ 
435 km). However, at these layers, in the intensity pattern only structures larger than 2 "0 are still connected with those at the 
continuum level although showing inverted brightness contrast. This confirms an inversion of temperature that we have found 
at a height of ~ 140 km. A possible evidence of gravity waves superimposed to the convective motions is derived from the 
phase shift analysis. We interprete the behaviour of the full width at half maximum and the equivalent width as a function of the 
distance to the granular borders, as a consequence of enhanced turbulence and/or strong velocity gradients in the inter granular 
lanes. 

> ' 

CO ■ Key words. Sun: general - Sun: photosphere - Sun: granulation - Techniques: spectroscopic - Hydrodynamics - Turbulence 
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1. Introduction First studies of the evolution of granulation were mainly 

connected with the determination of lifetimes. However, a 

O An important part of our knowledge about solar granulation detailed description of the evolution of individual granules 

CN ] has been obtained from the analysis of individual images demands data with much better quality than that required for 

■ and/or spectrograms. Actually, the most detailed images mere identification . Mehltretter (1978) using high resolution 
are obtained by applying different on-line and post-facto timeser i eso f w hit e light images found two types of behaviour: 
techniques (see e.g. Bonet 1999 for a summary). Additionally, Qn dme scales shorter than six minutes granules Qr groups of 

J ! one-dimensional (ID) and two-dimensional (2D) high spatial granules seem tQ expand At larger time scales the individual 

<S resolution spectrograms (« . 5) are used to investigate the granules fragment or drift Dialetis et aL ( 1986) describe d the 

dynamical behaviour of granules and intergranular lanes different types of formation (evo i ution of a sma n granu i e; 

and their penetration through the different layers of the fragmentation of a larger granule; and merging of smaller 

photosphere. A few parameters like continuum intensity, granules) and disa pp eara nce (fading into the background, 

line core intensity, equivalent width and full width at half merging a neighbouring granule or fragme ntation) of these 

maximum are often used, in addition to line bisectors, in order structures H irzberger et al. (1999) found that fragmentation 

to characterise the shape of photospheric lines and to infer and merging are the mQSt frequent mechanisms for birth and 

important physical parameters. Vertical granule velocities are disappearance respectively 
derived from Doppler shift measurements in the line core. 
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Concerning the spectroscopic methods, the vertical struc- 
ture of the convective motions has been studied in detail using 
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correlation and spectral analyses techniques applied to the fluc- 
tuations of different line parameters. For such studies different 
sets of observational data have been used. 

a) Most investigations have been based on single ID spectro- 
grams, Nesis et al. (1988, 1992, 1993), Komm et al. (1990, 
1991a, 1991b, 1991c), Hanslmeieret al. (1991a, 1991b, 1993, 

1994) , but without filtering properly the p-modes from their 
data, which strongly masks the convective signature. 

b) Tracking of spectra: Altrock et al. (1984) and Nesis et 
al. (1999, 2001) followed the evolution of selected granules 
during short periods of time. Keil (1984) and Johannesson 
et al. (1992) obtained 2D spectral images derived from ID 
spectrograms by scanning a certain area of the solar surface. 
More recently Collados et al. (1996) obtained improved data 
using the Correlation Tracker developed at the Instituto de 
Astrofisica de Canarias. 

c) 2D spectrograms have been obtained by using a 
Multichannel Subtractive Double Pass (MSDP) spectro- 
graph (Roudier et al. 1991; Espagnet et al. 1993, 1995) or 
Fabry-Perot interferometers (Salucci et al. 1994; Bendlin & 
Volkmer 1993; Krieg et al. 2000; Hirzberger et al. 2001). 

From all these investigations the granulation phenomenon 
can be currently described as follows: convective overshoot 
of the plasma flow from the solar convection zone into the 
photosphere forms a pattern of bright cellular elements show- 
ing upward motions - the granules, surrounded by a network 
of dark intergranular lanes, where downflow motions are 
observed. According to correlation analyses (Deubner 1988; 
Salucci et al. 1994; Espagnet et al. 1995) the horizontal tem- 
perature fluctuations associated with these motions decrease 
rapidly with increasing height until they vanish. Only the 
largest granules (> 1 "5) contribute to the brightness pattern 
observed above where the brightness contrast is inverted. Note 
that the values about the height where the intensity fluctuations 
vanish varies a lot in the literature. One finds values from 60 
km reported by Kneer et al. (1980), 60-90 km (Espagnet et al. 

1995) 170 km (Hanslmeier et al. 1993; Komm et al. 1991a) 
up to 270 km where Bendlin & Volkmer (1993) detected 
brightness signatures of the granulation. The values may differ 
from each other due to differences in the method to establish 
the geometrical height scale (transformation from T to z) 
and in different methods to filter the p-modes. The vertical 
velocity field persists in the upper layers, but waves (acoustic 
and gravity) and turbulent motions may also contribute to the 
velocity fluctuations (Deubner 1988; Deubner & Fleck 1989; 
Salucci et al. 1994; Espagnet et al. 1995). An overview on 
basic properties of solar granulation can be found for instance 
in Bray et al. (1984), Spruit et al. (1990), Muller (1999) and 
Brandt (2000). 

In order to get a detailed insight into the spatial and 
temporal evolution of granulation, a series of high spatial 
resolution spectroscopic observations is required, which is, 
however, extremely difficult to obtain. Some conditions should 
be fulfilled a priori for this study: 

a) A time series of spectrograms, long enough to allow an 
adequate filtering of the p-modes, because the presence of 



oscillations tends to decorrelate intensity and velocity fluctua- 
tions. 

b) In principle, a 2D coverage of the solar surface is desirable 
to sample the full evolution of granular areas. Nevertheless, 
time series of ID slit spectra still supply worthwhile infor- 
mation, provided we have an adequate tracking of the solar 
structures. 

c) A good seeing stability over the whole time series to allow 
the dynamical study of individual granules. 

We analyse a time series of ID slit spectrograms, recorded 
at the German Vacuum Tower Telescope (Observatorio del 
Teide, Tenerife) which fulfills the above mentioned conditions. 
In this paper we present a study of the global and spatial vari- 
ations of different line parameters in a quiet region at the solar 
disk centre. In Sect. 2 of this paper we describe the observations 
and the process of data reduction. In Sect. 3 the tools for data 
analyses are discussed, i.e. determination of line parameters 
and formation height of spectral lines, sub- sonic filtering, co- 
herence and phase spectral analysis, and image segmentation. 
In Sect. 4 global correlation between several line parameters 
derived from our spectrograms, for filtered and unfiltered data, 
are investigated. In addition, the mean values of line parameters 
calculated in granular and intergranular areas are compared. In 
order to get a deeper knowledge about the penetration of dif- 
ferent structures into the photosphere we perform in Sect. 5 a 
spatial coherence and phase shift analysis between the several 
line parameters and we study the dependence of these parame- 
ters on the spatial distance to the granular borders. In Sect. 6 we 
examine two particular granular cells, large and small, respec- 
tively, and the results are compared with the previous statistical 
findings for the whole sample. Finally, we summarise in Sect. 
7 the conclusions of this work. 

2. Observations and data reduction 

2.1. Observations 

Time series of CCD spectrograms were obtained in July 8, 
1993 with the 70 cm Vacuum Tower Telescope (VTT, see 
Schroter et al. 1985, for a description) at the Observatorio 
del Teide (Tenerife). From the recorded material, a time series 
spanning 50 min has been selected for this study, due to the 
stability of the seeing conditions. Figure 1 shows the time vari- 
ation of the continuum AI^ (normalised rms contrast of the 
granulation) along the slit. The strong spikes of the unfiltered 
signal are related with variations of seeing which also produce 
very high time-frequency structures (vertical narrow streaks) 
in the background image of Fig. 1. This information occurs in 
the k-co diagram within the domain removed by the p-modes 
filter. Consequently, the narrow streaks are smoothed out as a 
by-product of the p-modes filtering. This smoothing also af- 
fects the A/rms in the spatial direction (see dashed line in Fig. 
1). 

Magnetically active regions have been avoided by monitor- 
ing the entrance slit position in Call K slit-jaw images. The 
observations have been performed at the solar disk centre in 
order to obtain the vertical velocity fluctuations. 
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Line 


Element 


A (A) 


W X (mA) 


Pexc (eV) 


log gf 


Reference 


(I) 


Fel 


6494.994 


165.0 


2.39 




Moore et al. 1966 


(I) 


Fel 


6494.994 


162.0 


2.40 


-1.23 


Gadun et al. 1997, 1999 


(III) 


Fel 


6496.476 


69.0 


4.77 




Moore et al. 1966 


(III) 


Fel 


6496.472 


67.7 


4.79 


-0.66 


Gadun et al. 1997, 1999 



Table 1. Characteristic parameters of the lines analysed. The columns in order: line, element, wavelength, equivalent width, 
excitation potential, gf (g = statistical weight, / = oscillator strength) and reference to literature. 
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Fig. 1. Time variation of the raw AI^s along the slit (solid line) 
and after filtering of the p-modes (dashed line). Background: 
image of unfiltered continuum intensity in temporal (horizon- 
tal axis) and spatial direction (vertical axis). The total spatial 
direction covers a range of 46 "6. 



We have used a nitrogen cooled CCD camera with 1024 
x 1024 pixels. The exposure time was 1.5 s and the interval 
between two consecutive spectrograms 20 s. The entrance slit 
of the spectrograph was set to 100 ji (0 "46) and the CCD 
images included three photospheric Fel lines (A 6494.994, 
6495.740, 6496.476 A) that had been studied in previous 
papers (Hanslmeier et al. 1990) and will be hereafter referred 
to as I, II, and III, respectively (see Table 1 and Fig. 2). Line 
II has been finally excluded from this analysis due to a blend 
with a terrestrial water vapour line. A short summary of the 
most important properties of the spectral lines studied is given 
in Table 1 . The pixel size in the focal plane of the spectrograph 
was 0"093 in the spatial direction and 2.06 mA (« 95.15 
ms _1 ) in the spectral direction. In spatial direction, the slit 
covered a total length of 95 "2 (1024 pix) although after 
flatfielding, only the central part between the two reference 
hairs for positioning (46 "6 equivalent to 501 pix) was kept 
for further analyses. From the computed power spectra we 
estimate the effective spatial resolution achieved as ~ '.'5. 

To follow the time evolution of a physical point on the solar 
surface one must ensure a perfect tracking of this point on the 
spectrograph slit all the time. Image motion induced by atmo- 
spheric turbulence and mechanical vibrations of the telescope 
degrade the angular resolution and contaminate the observa- 
tion with information stemming from neighbouring points. 



Similarly, a drift of the image caused by non-perfect guiding of 
the telescope leads to a loss of tracking of the point of interest. 
To minimise these drawbacks originating from global image 
motion, our observations were performed through a Correlation 
Tracker developed at the Instituto de Astrofisica de Canarias 
(Bonet et al. 1992; Ballesteros et al. 1994; Bonet et al. 1994; 
Ballesteros et al. 1996) that was available for the first time in 
1993 for scientific observations. This device, mounted in front 
of the spectrograph slit at the VTT, has a tilting mirror which 
compensated for overall wavefront tilts - responsible of the 
global motion - at a fast rate. The system allowed a precise 
stabilisation of a portion of the solar image as well as the track- 
ing of the corresponding area on the solar surface for long peri- 
ods, following the solar rotation and the motions induced by the 
large-scale velocity pattern. This Correlation Tracker proved to 
have a good performance. Small faults in the tracking of the 
observed area, along the spectrograph slit (±1 px) were cor- 
rected a posteriori by cross-correlating the continuum intensity 
profiles of adjacent spectrograms in the spatial direction. 

2.2. Flatfielding 

The flat fields have been obtained in the same spectral region 
by moving the pointing of the telescope across the centre of the 
solar disk. The presence of sharp lines in wavelength direction 
(hereafter scratches) in the obtained spectrograms (see Fig. 2, 
left panel, where a particular raw spectrum is presented) and 
flat fields has required a special procedure to reduce our data, 
since parts of our further analyses would have been strongly 
affected by these artifacts. 

We have two types of scratches. The first stems from hair 
lines placed in front of the spectrograph slit (corresponding to 
the two darkest horizontal lines visible in the upper and the 
lower part of the raw spectrum and flat field source frame, re- 
spectively). The second type of scratches results from dust par- 
ticles on the spectrograph slit. 

In the following, we present a brief description of the data 
reductions we have applied (for a similar procedure see also 
Wohl et al. 2002): 

a) Subtraction of the dark current from the particular spectra 
and flat fields. 

b) Construction of the flat field source frame by averaging the 
individual flats, taken at the beginning and the end of the 
observing run. 

c) The removal of the spectral lines from the flat field source 
frame requires the division of the flat field source frame 
by its mean spectral profile averaged in spatial direction. 
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Fig. 2. Particular raw spectrum (left panel) and corresponding 
reduced spectrum (right panel). In both panels we can see from 
left to right: Line I, Line II, the telluric blend, and Line III. 



To perform this correction in a proper manner one has to 
rectify the possible inclination of the spectral lines with re- 
spect to the vertical edges of the frame. 

d) Splitting up the flat field matrix into two components: slit 
flat (containing the systematic effect induced by dust and 
wires on the spectrograph slit) and camera flat (including 
the rest of the spurious systematic effects) 

e) After dividing each spectrogram by the camera flat, we 
have compared by correlation the j-positions of the most 
prominent horizontal scratches in the slit flat and in each 
individual spectrogram. From this comparison the optimum 
shift of the slit flat to match the position of the scratches in 
the spectrogram has been derived with subpixel precision. 
Finally, each spectrogram has been divided by the displaced 
slit flat. A particular reduced spectrum is presented in the 
right panel of Fig. 2. 

f) Due to some remaining residuals of the hair lines in the re- 
duced spectrum, we have removed from each spectrogram 
the upper and lower parts outwards the hair lines, so that 
only the central area in the spatial direction is kept for fur- 
ther studies. 



2.3. Filtering of noise 

The noise and residuals of some remaining weak horizontal 
scratches in our spectra have been filtered out in the Fourier 
domain by applying a 2D optimum filter that is a generalisa- 
tion of the ID filter described by Brault & White (1971) (see 
also Bonet 1999). The mathematical expression for this filter 
(*(*)) is: 



*(s): 



ft(s)+/Ms)' 



(1) 



where s is the frequency, Ps(s) and P^(s) are smoothed models 
of the signal and noise power spectra averaged over all spec- 
trograms. Pn(s) has been modelled as a constant value equal to 
the mean level of the power in the range of high frequencies. 
The model for P$ (s) has been obtained by smoothing the mean 
power spectrum using a boxcar of empirically determined di- 
mensions and then subtracting P^(s). A 3D representation of 
the resulting filter is shown in Fig. 3. 



Fig. 3. 3D representation of the optimum filter to improve the 
signal to noise ratio and to eliminate the residuals of horizontal 
scratches in the spectra. 



Figure 4 shows a particular spectrogram before and after 
filtering. 

Each spectrogram has been also corrected for an incli- 
nation of the spectral lines with respect to the vertical edge 
of the CCD. Finally, the average continuum of each spectro- 
gram has been fitted to the local continuum of the Liege Atlas 
(Delbouille et al. 1973-1979) at this wavelength range. 

3. Techniques of data analysis 

3. 1. Computation of line parameters 

Line parameters have been calculated for Line I and III from 
each spectrogram of our time series as follows: 

The continuum intensity (hereafter I con ) has been deter- 
mined as the average intensity value of 5 points, placed around 
the position of the maximum value of intensity in the right 
wings of Line I and Line III, respectively. 

The position of the minimum of the lines has been deter- 
mined by fitting a fourth order polynomial to 12 points (Line 
I), 10 points (Line III) and 8 points (terrestrial line), at the bot- 
tom of each line. By taking the intensity value at this position, 
the line core intensity of the solar lines (hereafter I\ and Tni) is 
obtained. 

To measure the line core velocities the following procedure 
has been used. All spectra have been averaged in space and 
time, thus obtaining a mean profile on which the minimum po- 
sitions of the two solar and the terrestrial lines have been com- 
puted. The absolute distance in wavelength between Line I and 
the telluric line, d\, and between Line III and the telluric line, 
dm, has been adopted as the reference zero- velocity distance. 
This assumption is justified because we are only interested in 
relative velocity variations and not in absolute velocity values. 
The evaluation of the minimum positions has been performed 
also on each individual profile. By calculating the distances 
from Line I and Line III to the telluric line, d\ and dm, the 
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Reduced Spectrum Filtered Spectrum 




Fig. 4. 3D representation of an individual spectrogram before (left panel) and after (right panel) applying the optimum filter to 
improve the signal to noise ratio. Spectral-axis: wavelength (mA) in relative units; spatial-axis: arcsec\ vertical-axis: intensity in 
relative units. 



Doppler velocities of each line have been derived as the differ- 
ences d\ — d\ and dm — dm- Instead of converting the Doppler 
shifts into velocities hereafter the velocities will be referred in 
terms of distances, that is, V\ = d\ — d\ and Vm = dm — dm. 

Because of a blend in the upper left wing of Line III, the 
equivalent widths of both solar lines (hereafter EW\ and EWm) 
have been computed by integration of the area enclosed by each 
line profile below the 90%-level of the continuum. 

Additionally, the full widths at half maximum for Line I 
and III (hereafter H W\ and HWm) have been computed. 

In the following we will refer to /, V, EW, HW, (without 
subscripts) when a common aspect in both lines is described. 
The values of the different line parameters (i.e. continuum in- 
tensity (7 C on), line core intensity (/), line core velocity (V), 
equivalent width (EW) and full width at half maximum (HW)) 
at each point x along the spatial direction are arranged vs. time 
(t) in a 2D representation that from now on will be termed "im- 
age" (e.g. image of intensities, equivalent widths, etc.). 

3.2. Filtering ofp-modes 

The influence of p-modes in our images has been minimised 
by applying a subsonic filter in the Fourier space (k x , CO). Only 
those spectral components below the straight line CO = v p k x are 
retained whereas the rest are set to zero (Deubner 1988). The 
constant v p defines the cutoff, i.e. the maximum phase velocity 
admitted by the filter. Values of v p ranging from 3 km s _1 to 6 
kins" 1 have been tested obtaining similar results, and v p = 5 
km s _1 has been finally taken as a conservative choice. 

Since the granular brightness fluctuation is a 2D spatial pat- 
tern, a 3D description, g(x,y,t), would be desirable to perform 
a more complete spectral analysis -here g means continuum 
intensity or any other spectroscopic parameter. However, our 
spatial information is restricted to only one dimension (along 
the spectrograph slit, x-axis, perpendicular to the dispersion di- 
rection), so that we have been compelled to work in 2D, (x,0> 



and consequently to obtain the Fourier transform G(k x ,co) 
which represents the integration of the 3D Fourier transform, 
G(k x ,k y , co), over k y . 

In effect, let g(x,y,t) be the 3D array consisting of a time 
series of 2D images describing granulation parameters (inten- 
sities, Doppler velocities, equivalent widths, etc). The relation 
between G and G is derived from the inverse Fourier transform 
formula for the case y = 0, i.e. g(x,0,t): 

g(x,0,t) = J fG{k Xl co)e 2Ki ^ m Uk x dco, (2) 

where 

G(k x ,co) = J G(k x ,ky,co)dk y . (3) 

According to Eq. 3, the information at a given k x is con- 
taminated by the contribution from spectral components lying 
onto a parallel to the k y -axis at this specific k x . This merging of 
spectral information evidences that subsonic filtering would be 
desirable in G(k x ,k y , co) instead of in G(k x , (o), the filter being 

defined by co = v p k, where k = ^jkl + k] (see Title et al. 1989). 

For testing to what extent the 2D treatment is a reliable 
approximation in the filtering process we have performed a 
numerical experiment by using an 80 min time series of 330 
x 330 pixel high resolution images of granulation taken at 
the SVST (Swedish Vacuum Solar Telescope, Observatorio 
Roque de los Muchachos, La Palma), with a pixel size of 
"062 x "062 and a time spacing of 18.9 s, see Hirzberger et 
al. (1997). A 2D array (x, t) has been extracted from the original 
3D array by fixing the y-coordinate in all granulation images; 
the resulting array emulated the one obtained for the contin- 
uum from the VTT's time series of spectrograms (that we have 
termed "image" of continuum intensity). We also have aver- 
aged over 8 pixel in the y-direction to simulate the slit width 
of the spectrograph (cf. Sect. 2.1). Then, 3D and 2D subsonic 
filtering with a common cutoff have been applied to the orig- 
inal 3D and to the extracted 2D array, respectively. We have 
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repeated the 2D extraction at the same v-position but this time 
from the filtered 3D array, and compared with the results of 
the 2D filtering. The match of the results has improved when a 
gradual rather than an abrupt cutoff has been considered in the 
2D filter. A cosine bell above the line CO = v p k x , over the 10% 
of the length of the co-coordinate has been applied. Figure 5 
shows a representative example of the results for the case of 
one scan in the x-direction. Note that the differences between 
the curves from 3D and 2D filtering - thin and thick solid lines, 
respectively - are small as compared with the original unfil- 
tered curve - dotted line. 
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Fig. 5. Intensity profiles of solar granulation along the spatial 
coordinate x: Original (dotted) and after 2D (thick solid) and 
3D (thin solid) subsonic filtering. 



3.3. Calculation of coherence and phase difference 
spectra 

Analysis of coherence and phase differences between intensity 
(0 and velocity (v) fluctuations of different spectral lines pro- 
vides an adequate diagnostic to discriminate between the var- 
ious velocity fields that contribute to the observed dynamics 
in the solar photosphere. Thus, [v-v], [i-i] and [v-i] coher- 
ence and phase difference spectra are usually calculated with 
respect, separately, to both wavenumber, k x , and temporal fre- 
quency, co (see Deubner & Fleck 1989 and references therein). 

Let g\(x,t) and g2 (x, t ) be 2D arrays representing the varia- 
tion in space and time of two granulation parameters. We have 
calculated the coherence, p\2, and the phase difference, 0i2, 
spectra between these two arrays, with respect to k x , as follows 
(see Edmonds & Webb 1972; Deubner 1988) 



Pl2(k x ) 



£Gi(fe,0G5(fe,0 



£ Gi(M)GKM) 



,Im£Gi(M)Gl(M)> 

^12 fe) = arctan — =^ 

121 x) ^Re£Gi(fe,0(^(fe,0 



(4) 



(5) 



where * stands for complex conjugate, G\ and G2 are the 
Fourier transforms of g\ and g2 with respect to the spatial vari- 
able and the summations are taken over times t . An interchange 
in the role of the spatial and temporal variables allows to obtain, 
similarly, coherence and phase difference spectra with respect 
to CO. 

3.4. Heights of formation of the spectral lines 

We aim at studying the variation of some physical parameters 
with height in the photosphere. Therefore, it is appropriate to 
determine how and where the two spectral lines respond to 
the temperature and velocity field fluctuations. The response 
functions (RF) describe the sensitivity of the intensity profile to 
different photospheric quantities, like for instance temperature 
T or line of sight velocity Vlos (Ruiz Cobo & del Toro Iniesta 
1994; Sanchez Almeida et al. 1996). We have calculated 
the RF of Line I and III to T perturbations (Rt) and Vlos 
perturbations CRv L0S )> applying the SIR code (Ruiz Cobo & del 
Toro Iniesta 1992). We refer also to Sheminova (1998), Gadun 
et al. (1997, 1999) and Perez Rodriguez & Kneer (2002) for 
other calculations of response functions. 

Line I is strong enough to consider that NLTE effects may 
cause significant errors in the evaluation of the RF. We have 
used departure coefficients j3i ow and j3 up obtained by Shchukina 
& Trujillo Bueno for the quiet sun model of Maltby et al. (1986, 
MACKKL model) and a 63 levels model atom of neutral iron. 
Aow (ftp) stands for the ratio between the population of the 
lower (upper) atomic level evaluated in NLTE and LTE con- 
ditions, respectively. The absorption coefficient ^ NLTE and the 
source function S^j LTE can be written as: 



K 



NLTE 




ftp/ftow<? 



-hv/kT 



-hv/kT 



(6) 



nNLTE 



2hV 5 



1 



C L 

2hv 3 



Aow//V w - 



-hv/kT 



fto 



ftp gLTE 



ft 



(7) 



The approximation in the second term of these equations 
is valid when we are able to neglect the stimulated emission 
(hv » kT), like in the case of solar lines in the visible range 
of the spectrum. In Fig. 6, j3i ow , /3 up and its ratio is plotted ver- 
sus T, the continuum optical depth at 5000 A. Using Eq. 6, it 
is clear that the absorption coefficient is decreased by NLTE 
effects in the range of log T = to — 4 and increased from log T 
= —4 on, respectively, because of the behaviour of j8i ow . The 
NLTE source function shows a strong decrease with respect 
to the Planck function from log T = — 3 on caused by the be- 
haviour of the ratio j3 up / j3i ow (see Eq. 7). 

By substituting ^ NLTE and S^ LTE for K and S y in the 
expression of the response function (see Ruiz Cobo & del Toro 
Iniesta 1994) we have an suitable approximation of the RF 
which neglects changes in the departure coefficient but takes 
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departure coefficients 



10.0 =- 




Fig. 6. Departure coefficients of Line I calculated for the quiet 
sun model of Maltby et al. (1986, MACKKL model) and a 
Fe -atomic model with 63 levels, where j3i ow and /3 up are the 
departure coefficients of the lower and upper atomic level, re- 
spectively. The values have been kindly provided by Natasha 
Shchukina & Javier Trujillo Bueno. 




into account the more important NLTE effects (for a discussion 
see Socas-Navarro et al. 1998). 

To evaluate the RFs we have used the MACKKL model, 
assuming a constant macroturbulence £ mac = 1.7 kms -1 , a 
constant microturbulence ^mic = 0.9 kms -1 , a collisional 
damping enhancement factor E = 1.3, and a loggf value = 
— 1.23 for Line I and values of £ mac = 1.3 kms -1 , 4^ = 0.9 
kms -1 , E = 2.9 and log gf = —0.66 for Line III, respectively. 
With these selection of parameters the synthesised profiles 
match the mean observed ones. For the calculation of the 
geometrical heights scale of this model we have solved the 
hydrostatic equilibrium equation considering a boundary con- 
dition for the surface gas pressure of 0.36 dyn/cm 2 (MACKKL 
model). 

In the left panels of Fig. 7 we can see the RFs to T 
perturbations for Line I (top) and Line III (bottom) and in 
the right panels the corresponding response functions to Vlos 
perturbations, respectively. The RFs plotted in this figure stand 
for the change of intensity at each wavelength when a local 
perturbation of 100 K and 1 kms -1 (the order of magnitude 
of T and Vlos perturbation we expect in granulation) is intro- 
duced at each layer in log T. The two selected spectral lines are 
well differentiated in altitude and span a large height/optical 
depth range of the photosphere. 

Following Ruiz Cobo & del Toro Iniesta (1994) the RFs of 
the line parameters /, V, EW and H W to T and Vlos perturba- 
tions can be straight forwardly expressed in terms of the RF of 
intensity. In Fig. 8 we have plotted the RF of /, EW, HW for T 
perturbations and the RF of V for Vlos perturbations. Each line 
parameter is particularly sensitive to one photospheric quan- 
tity. However, all the different line parameters are sensitive to 
photospheric quantities in a broad range of the photosphere. 



Fig. 7. Plot of the RF for T (left panels) and Vlos perturbations 
(right panels) for Line I (upper panels) and III (lower panels), 
where Rt and ^y L0S are measuring the change of intensity for T 
perturbation of 100 Kelvin and Vlos perturbations of 1 km s -1 
at each layer of width A log t = 0.1. 
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Fig. 8. Panels a, b and c correspond to the RF of I, EW and 
HW to a T perturbation of 100 K. Panel d shows the RF of V 
to Vlos perturbation of 1km s -1 at each layer of width A log T 
= 0.1. RFs of Line I and III are represented as solid and dashed 
lines, respectively. 



Especially difficult is the case of the RF of V\, which shows 
a maximum sensitivity in the range spanning from log T c± — 1 
to log t ~ —3.1, so that we are compelled to work with sam- 
pled information originating from very different heights of the 
photosphere. 
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Line [log t] / EW HW V_ 



I 


range 


[0., -5.] 


[l.,-5.] 


[0.5, -5.] 


[0., -5.] 




H0F 


-3.1 


-0.6 


-0.7 


[-1.2, -3.1] 


III 


range 


[0., -3.] 


[1., -3.] 


[1..-3.] 


[0., -3.] 




HOF 


-0.9 


-0.6 


-0.4 


-0.8 



Table 2. Column 3, 4, 5, 6: total range of sensitivity to T and 
Vlos perturbations and height of formation (HOF) of /, EW, 
HW and V in optical depth, for Line I and III, respectively. 



Line 


[km] 


/ 


EW 


HW 


V 


I 


range 


[0, 1065] 


[-66, 1065] 


[-36, 1065] 


[0, 1065] 




HOF 


435 


75 


90 


[166, 435] 


III 


range 


[0, 421] 


[-66, 421] 


[-66, 421] 


[0, 421] 




HOF 


120 


75 


46 


105 



Table 3. Column 3, 4, 5, 6: total range of sensitivity to T and 
Vlos perturbations and height of formation (HOF) of /, EW, 
HW and V in geometrical height, for Line I and III, respec- 
tively. 



The RF of EW is mainly influenced by the effect that in 
higher layers (above log T = 0) an increasing temperature yields 
a negative lobe in the RF caused by a smaller absorption due to 
the decrease of Fe I abundance (Saha formula), while in lower 
layers (below log % = 0) a positive lobe in the RF is caused by 
an increase of continuum intensity. 

The RF of HW is influenced by the same effect like in the 
case of EW plus an additional broadening caused by the in- 
crease of Doppler width and resulting in a second positive lobe 
in the response function in higher layers. Note, that the height 
of formation for EW and HW would change, if we would use 
other definitions like the centre of gravity of the response func- 
tion instead of its peak. 

The sensitivity of /, EW and HW to Vlos perturbations 
and of V to T perturbations is negligible, because our model 
did not represent strong velocity gradients. 

In Tables 2 and 3 we have summarised the results regard- 
ing the height of formation. Total ranges of the sensitivity 
of different line parameters to T or Vlos perturbations and 
assumed heights of formation are given in optical depth and 
geometrical height. 

The formation height concept is somewhat inexact and to 
ascribe certain physical properties derived from spectral lines 
to a given height in the atmosphere is maybe misleading. 
However, it is useful for interpreting the results qualitatively. In 
subsequent papers the observed data will be inverted by means 
of using the complete response functions to the different physi- 
cal quantities. This will allow us to obtain results free from the 
simplification assumed here. 

3.5. Selection of granular and intergranular regions. 

We must first consider how to separate granules from inter- 
granular lanes. In general we can formulate the problem in this 
way: granules (intergranular lanes) are defined by all those 
pixels (i, j) with intensity and velocity values satisfying the 



following condition: 
Granules: 

7ij> (l+q)I ; V ij >(l+q)V (8) 
Intergranular lanes: 

4j< (l-q)7 ; Vii<{l-q)V (9) 

where /y and Vy are intensity and velocity values at each 
pixel in space and time, and 7, V the mean values of brightness 
and velocity over the whole image, respectively, q and q are 
threshold parameters defining a granular-intergranular transi- 
tion region. 

Most of the authors have only used an intensity criterion 
by assuming q = 0. Keil (1984) defined granules (intergranular 
lanes) as those regions where the difference between the local 
continuum intensity and the mean continuum intensity (7) is 
above (below) 8 % (q = 0.08). This percentage was reduced 
to 1% by Muller & Roudier (1984) and to 4 % by Karpinsky 
(1989). 

The granulation as a convective phenomenon is better de- 
tected in the deepest observable layers of the photosphere. 
The continuum intensity reflects the temperature fluctuations in 
these layers, whereas line core velocities of Line III and Line 
I stem from higher layers. Since we do not have information 
about velocities from the deepest layers it seems more reason- 
able to use the single criterion based on intensity. Thus, finally 
this single criterion has been adopted, assuming q=0. 

4. Global correlations 

4.1. A glance to the images 

Figure 9, composed of nine vertical panels, summarises some 
of the results of our VTT spectrogram series. Each image, com- 
posed of 152 columns, represents the variation of a spectro- 
scopic parameter along the spectrograph slit (vertical), at suc- 
cessive times throughout the series (horizontal). The contents 
of the individual panels or images are as follows: 

• 1 : binary image of the continuum brightness showing the 
time evolution of the granular (white) and intergranular 
(black) regions along the spectrograph slit. 

• 2 and 3: unfiltered and filtered continuum brightness (/con)- 
The very fine vertical structures (high temporal frequency) 
observed in our raw images, caused by differential seeing 
distortion between successive spectrograms or by a residual 
tracking error are removed by the subsonic filter. 

• 4 and 5: unfiltered and filtered line core velocities for Line 
HI (Vm). 

• 6 and 7: unfiltered and filtered line core velocities for Line 
I (Vi). 

• 8 and 9: filtered line core intensities for Line III and I (/m, 
/i), respectively. 

From these images we can extract the following qualitative 
information: 
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3 4 5 6 7 8 9 




25 50 
time [min] 

Fig. 9. Images representing the variation of continuum brightness, velocities and line core intensities for Line I and Line III. From 
the left to the right: (1) 7 con binary image, (2) 7 con unfiltered, (3) 7 con filtered from p-modes, (4) Vm unfiltered, (5) Vqi filtered, 
(6) V\ unfiltered, (7) V\ filtered, (8) filtered, (9) I\ filtered. Horizontal axis is time and vertical axis space. To facilitate a visual 
comparison between velocity and intensity images, we have reversed the scale of velocity values so that granular velocities 
(negative) are represented also as bright areas. See the text for a detailed description of the different panels. 





Vm 


Vi 


/m 


h 


EWm 


EWi 


HW m 


HWi 


/con 


unfiltered 


-0.61 


-0.17 


0.14 


-0.22 


0.50 


0.34 


-0.43 


-0.37 




filtered 


-0.87 


-0.49 


0.03 


-0.44 


0.69 


0.58 


-0.46 


-0.37 


Vm 


unfiltered 
filtered 




0.74 
0.62 


0.02 
0.04 




-0.48 
-0.79 




0.25 
0.30 




Vi 


unfiltered 
filtered 








0.07 
0.48 




-0.18 
-0.51 




0.02 
0.24 


/m 


unfiltered 
filtered 








0.32 
0.24 


-0.54 
-0.44 




0.63 
0.68 




k 


unfiltered 
filtered 












-0.83 
-0.78 




0.56 
0.56 


EW m 


unfiltered 
filtered 












0.79 
0.82 


-0.31 
-0.33 




EWi 


unfiltered 
filtered 
















-0.09 
-0.03 


HW m 


unfiltered 
filtered 
















0.45 
0.54 



Table 4. Correlation coefficients between the fluctuations of different line parameters for filtered and unfiltered data. 



• The filtering method has apparently removed the fluctua- 
tions associated with the p-modes. 

• The expected difference between granules and intergranu- 
lar lanes can be clearly seen. The width of the structures 
is also the one expected (cf. Bray, Loughhead & Durrant 
1984). There exist granules and intergranular lanes that 
persist over the full observing period in agreement with 
Sobotka et al. (1999). 

• In the different images we see small-scale structures that 
could correspond to fine features inside the granules or per- 
haps to non-convective motions present in the upper layers. 

• As derived from numerical experiments we have done, the 
small amplitude cross-shaped structure visible in images 7, 
8 and 9 is caused by the cutoff of the p-modes filter, even 
in the case of smoothing the transition from 1 to 0, within a 
narrow band close to the cutoff. 

• The granules suffer different kinds of evolutionary events 
like expansions, merging and fragmentation as it is ob- 
served in 2D images (e.g. Hirzberger et al. 1999). 



4.2. Tentative explanation of global correlations 

The filtering process allows us to discriminate between two 
different domains, namely: a) the oscillatory component, 
representing the p-modes together with residual motions 
induced by the tracking error and/or seeing, and b) the connec- 
tive component, representing the signature of the convective 
overshooting in the photosphere together with the gravity 
waves (Deubner & Fleck 1989; Bonet et al. 1991; Komm et 
al. 1991c). We have calculated the global correlations between 
different line parameters in order to get some preliminary hints 
about how the dynamics at different heights change through 
the photosphere and we have summarised these results in Table 
4. In most of the cases the filtering of p-modes in our images 
has increased the degree of correlation. In Sect. 5 we will 
enter in more detail by performing a coherence and phase shift 
spatial analysis to obtain information about the penetration of 
different structures across the photosphere. 

In the following we will remark the most significant results 
of correlation coefficients presented in Table 4 and give a ten- 
tative explanation in terms of the RF plotted in Fig. 8: 
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unfiltered filtered 




I 

con 



Fig. 10. Scatter diagrams of continuum intensity and velocity 
as measured from Line III (upper panels) and Line I (lower 
panels) before (left panels) and after filtering of p-modes (right 
panels). 



It also explains the weak positive correlation between the 
pair 

We do not find any significant correlation between the pairs 
t^con - ^iiiL Mil - ^111] and [Ii - /nil- This could be caused 
by the fact that the peak of the temperature response func- 
tion of Iiu is placed ~ log t = —0.9 (panel a of Fig. 8) just 
at the layer where the inversion of temperature occurs: for 
the case of a granule, higher temperatures in deeper layers 
compensate the variation of / introduced by lower temper- 
atures in the layers above. 

Berrilli et al. (1999), by studying two lines with formation 
heights about 68 km and 370 km, found correlations 
between line core velocity and line core intensity of 
opposite sign (positive for the line whose core is formed at 
~ 370 km). It is well known that the granular brightness 
fluctuations decrease rapidly with increasing height, vanish 
at ~140 km where the inversion of temperature contrast 
takes place, and increase slightly above this height again 
(Rodriguez Hidalgo et al. 1996). If we take into account 
that Iiu and I\ can be ascribed to a height of ~120 km 
and ~435 km (see Table 3), respectively, our results agree 
with the results obtained by Berrilli et al. (1999). Note, 
that our results are also in agreement with other works 
in the literature but there exist discrepancies concerning 
the height where the brightness fluctuations vanish (see 
Introduction). 



• Correlations with velocities: 

In Fig. 10 we present a scatter plot of / con vs. V, before and 
after filtering the p-modes. The high degree of correlation 
between 7 con and Vqi after filtering the p-modes is a clear 
consequence of the convective nature of granulation. The 
correlation between 7 con and V\ is weaker , thus reflecting 
a less efficient penetration of the convective velocity field 
into higher layers (Deubner 1988; Espagnet et al. 1995). 
This behaviour is also confirmed by a decrease of the coher- 
ence between [7 con — Vj] with respect to that of [7 con — Vm] 
(see Fig. 11). Remember, that from the calculation of re- 
sponse functions in Sect. 3, the height of formation (HOF) 
of Vm can be assigned to log T = —0.8, while Vi originates 
from a broad range in optical depth [—1.2; —3.1], clearly 
higher in any case. Thus we are able to conclude that con- 
vective motions penetrate up to the highest layers of the 
photosphere, where the velocity pattern is still connected 
with the brightness pattern at the continuum layer. 
The correlation between velocities at the two levels 
(Vm vs. Vi) is high but decreases considerably when 
we study only the convective component. This means, 
that the contribution of oscillatory velocities due to the p- 
modes is significant at the height of formation of both lines. 

• Correlations with line core intensity: 

We find a weak negative correlation between the pair [7 con 
- Ii] due to the reversal of temperature contrast in high pho- 
tospheric layers (e.g. Deubner 1988; Collados et al. 1996; 
Rodriguez Hidalgo et al. 1996). The inversion of tempera- 
ture contrast was predicted by theoretical models (Steffen et 
al. 1989; Stein & Nordlund 1989; Gadun et al. 1997, 1999). 



• Correlations with equivalent width: 

For the pairs [7 con — EW\\ and [7 con — £Wqi] we find for 
the convective component significant positive correlation. 
From panel b of Fig. 8 we can deduce that granular line 
profiles have larger equivalent width than the intergranu- 
lar ones, because granules are much hotter below log T = 
(where the RFs for both lines are positive) and slightly 
cooler above log X— — 1 (negative RFs). Also Hanslmeier 
et al. (1993) reported a positive correlation (~ 0.53) be- 
tween these parameters studying four lines, having a line 
core formation height between 70 and 150 km, although 
their data were not filtered. The positive correlation has 
been explained (e.g. Gadun et al. 1997) as due to the in- 
fluence of oscillatory motions. Our data show the opposite 
behaviour since we find enhanced correlation after filtering . 
The negative correlation between the pairs [Vqi — £Wqi] 
and [Vi — EW\] can also be explained taking into account 
the behaviour of equivalent width as stated in the previous 
paragraph. 

For the pairs [Im — EWm] and [I\ — EW\] we have found for 
the convective component a high negative correlation for 
Line I and a weak negative correlation for Line III. This can 
be explained taking into account the behaviour of line core 
intensity that we comment in case of correlation between 
7 con and /. 

The equivalent widths of both lines are highly correlated. 
This can be explained by the fact that the temperature 
response functions for the two lines are very similar to 
each other (panel b of Fig. 8). An increase in temperature 
produce a similar change in equivalent width resulting in a 
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strong correlation for the pair [£W[n — EW{\. 

• Correlations with full width at half maximum: 

The RF of HW for Line I has a strong negative lobe around 
log T = — 1 . As temperature fluctuations at these layers do 
not affect the continuum intensity, no significant correlation 
between 7 con — HW can be expected. For Line III, a weak 
correlation between these parameters is observed, because 
the negative lobe of the corresponding response function 
is placed at deeper layers. Besides, a weak correlation for 
the pair [//Wqi — HW\] can be expected, as the response 
functions of both lines partially coincide. 
Another case is the observed correlation between line core 
intensity and half width. Turbulence and strong velocity 
gradients produce an enhancement of line width together 
with an increase of line core intensity, while a growing 
of line width induced by temperature changes always goes 
with a decrease of line core intensity. Thus, the significant 
correlation between / and HW indicates that either turbu- 
lence or strong velocity gradients must be important agents 
to explain the enhancement of the width of both lines. Nesis 
et al. (1999) related HW with an unresolved turbulent ve- 
locity field placed mainly at the granular borders. 

We have not found any significant values for the remaining cor- 
relations between the several line parameters. 

4.3. Line parameters in granular and intergranular 
regions 

Mean values of the relevant line parameters were calculated 
separately for granular and intergranular regions (Table 5). 
As expected from the analysis of the correlation between the 



Mean values 


Granular regions 


Intergranular Lanes 


Ii 


0.26 ±0.01 


0.27 ±0.02 


Ira 


0.50 ±0.02 


0.50 ±0.02 


EWi 


116.1 ±2.8 


112.7±3.5 


EW m 


41.3±1.2 


39.4±1.5 


HW] 


186.6±3.6 


189.2±4.6 


HWm 


116.8±4.4 


120.6±4.8 



Table 5. Mean filtered values (with standard deviations) of line 
parameters in granular and intergranular regions. Granular (in- 
tergranular) regions are defined as those where 7 con is above 
(below) its mean value. 



several line parameters, we do not find significant difference 
between the mean values of line core intensities above granules 
and intergranular lanes in case of Line III, since the informa- 
tion contained in Tni originates mainly from layers around the 
inversion of temperature where the brightness contrast between 
granules and intergranular lanes remains near 0. For Line I we 
find a slightly lower mean value of line core intensity above 
granules than above intergranular lanes, thus we can conclude 
that in higher layers granules appear slightly cooler than the 
surrounding intergranular lanes due to the inverted tempera- 
ture contrast. This is in agreement with models of Steffen et 



al. (1989) and Stein & Nordlund (1989), which predict that 
temperature fluctuations should change sign in the overshoot- 
ing layer, which should be cooler above the ascending granular 
parts and 'relatively' warmer above the descending intergranu- 
lar regions. 

For both lines we find enhanced values of EW above gran- 
ules. Gadun et al. (1997) predict an increase of EW in granules 
as observed. HW shows higher values above the intergranular 
lanes, which has been also reported by Nesis et al. (1990) and 
Hanslmeier et al.(1993). 

5. Spatial variations 

5. 1. Coherence and phase shifts spatial spectra 

An analysis of coherence and phase difference spectra between 
two photospheric lines formed at well separated heights gives 
information about how the granulation phenomenon is related 
at these two layers (e.g. penetration of granulation through the 
photosphere). In Fig. 11, the coherence and phase difference 
spectra between several line parameters are plotted. Only phase 
shifts with a corresponding coherence above ~ 0.5 will be con- 
sidered as reliable enough. Coherence and phase shift are repre- 
sented vs. spatial frequency and result from temporal averages. 
In the present discussion we are assuming the upper limit of the 
granular scales at w 4 "0. 

First we want to remark that the computed coherence spec- 
tra increase significantly in case of the filtered images due to the 
elimination of p-modes. This is especially noticeable in case of 
Line I, because the amplitude of p-modes increases with height 
in the photosphere. 

From Fig. 11 we can infer the following ranges of spatial 
scales with significant coherence (i.e. > 0.5): a) scales above 
~ 1 "0 for [/ con - Vi], [7 con - V m l and [Vi - V m ]; b) scales be- 
tween - 9 "0 and 2 "1 for [/ con - / T ]; and c) scales > 1 "6 for 
[Vi — Zt]. These high values of coherence reveal that the gran- 
ular phenomenon at larger scales penetrates up to the highest 
levels of the photosphere (at least up to ~ 435 km). Lets study 
more deeply these three cases. 

• [/con - ViL [/con - VmL and [Vi - V m \. 

7 CO n leads Vni by 1 80° and Vi by a decreasing amount from 
180° to ~ 150° in the granular range with high coherence 
values (1 "0 - 4 "0). For larger scales (> 4 "0) the phase 
difference between / con and V\ continues decreasing. 
Consequently with these results V\ leads Vm by an amount 
augmenting from 0° to ~ 20° in the same granular range. 
Deubner (1989) finds similar results by studying the lines 
Cl 5380 and Fel 5383, with corresponding heights of 
formation ~ and 200 km above T5000 = 1- Similarly to 
our case, he finds in the granular regime the Fe I velocity 
signal leading the Cl velocity by 10° to 15°; however at 
the larger scales his velocities are well in phase. 

• [/con-/l]: 

The coherence at very large scale is below 0.5, then rises 
above 0.7 at about 5 "0 and drops abruptly at ~ 2 '.'\. Thus, 
regarding the brightness structure only the largest granular 
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Fig. 11. Spatial variations : Average behaviour of the coher- 
ence (coh) and phase shifts (phd) between the pairs [7 con — Vm], 
[/con - ViL [Vi - V m ], [/ con - /fl and [Vi - hi before (dashed- 
dotted line) and after filtering (solid line) the p-modes. 



difference similar to that one that we have encountered be- 
tween Vm an d Vi as caused by the increase of thermal relax- 
ation time with height and also by the presence of gravity 
waves superposed on the convective motions. 

We have not found any coherence for the pairs [7 con — hn] , 
[Vqi — hn] and [I\ — Im] due to the small brightness contrast 
at layers close to the inversion of temperature, resulting in 
a decrease of the signal to noise ratio and a weakening of 
coherence. 

Our results complement those by Deubner (1989) in the 
sense that we report information at complementary heights in 
the solar photosphere. Espagnet et al. (1995) perform a phase 
and coherence spectral analysis between intensity and veloc- 
ity fluctuations as measured at several positions (wavelengths) 
throughout the profile of the NaD2 line (A 5890) which corre- 
spond to various heights in the photosphere ranging from the 
continuum level up to ~ 550 km. In contrast with our findings 
they conclude that temperature fluctuations of granulation do 
not penetrate higher than about 60-90 km. By comparing the 
intensity fluctuations at the continuum level with those at dif- 
ferent levels in the photosphere, they find high coherence for 
large granules (>1"4) up to 60 km, vanishing at 90 km. For 
Nesis et al. (1988) and Komm et al. (1990) the coherence of 
granules > 1 '.'A vanishes at 170 km above the continuum level. 

Regarding velocity features, Espagnet et al. (1995) com- 
pare the velocities at different layers with both, the continuum 
intensity and the velocity at the height of 30 km. From these 
coherence and phase shift analyses they conclude that the over- 
shooting velocities of granules > 1 "4 and > 1 "6, respectively, 
cross the whole photosphere up to a level of at least 550 km, 
and that the height of penetration decreases with decreasing 
size. These results are in a closer agreement with ours except 
for two aspects: 1) they find almost perfect antiphase at all pho- 
tospheric layers, whereas we obtain that 7 con leads V\ by a vary- 
ing angle between 180° and ~ 150° (this decreasing trend in 
the phase difference for large granules could be a signature of 
gravity waves in the upper photosphere); and 2) and velocities 
of granules larger than 1 "0 instead of 1 "4 penetrate into high 
layers (435 km). 



scales (<~ 2 ^0 - 4 "0) are detectable at highest photo- 
spheric levels. Nevertheless, the brightness pattern at these 
layers is inverted with respect to that of the continuum 
as a consequence of the temperature inversion at ~ 140 km. 

• [Vi-h]: 

At the highest levels the brightness pattern is similar to 
that of velocities but only at the largest granular scales (~ 
1 "6 - 4 "0). From the phase shifts measured for [7 con — 1\] 
and [/ con — Vi], we would expect V\ leading I\ by an in- 
creasing value from 0° at 1 '.'6 to ~ 20° at 4 "0. For Fe I 
5383, Deubner (1989) finds the line core intensity lead- 
ing the velocity in a similar amount, although with an in- 
creasing trend toward mesogranular scales which can also 
be detected in the very low frequency range of our figure. 
Deubner (1989) explains this effect as well as the phase 



5.2. Dependence of line parameters with the distance 
to the granular border 

In this section we intend to figure out how the line parameters 
vary as a function of the normalised distance to the granular 
border. For this purpose granules and intergranular lanes have 
been identified in the spatial direction (columns in the image 
of the continuum intensity) defining the granular borders at the 
inflection points of the intensity profiles. The middle position 
of every interval between two consecutive inflection points is 
taken as the centre of the corresponding granule or intergranu- 
lar lane. 

For each granule the distances from the border to the cen- 
tre have been computed, normalised to unity and interpolated 
to the same number of distance steps, so that, independently 
of the granular size, the distance is zero at the border and 1 
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at the centre and containing 100 values in between. Similarly, 
for intergranular lanes the distances have been set to vary from 
zero at the border to — 1 at the centre, independently also of the 
intergranular lane size. 

Figure 12 shows the averaged values of the different line 
parameters for each bin of 0.02 normalised distance. 
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Fig. 12. /con, V, I, EW and HW vs. normalised distance to the 
granular border for small and large granules (G) and intergran- 
ular lanes (IG) for the filtered images of both lines. Small struc- 
tures (0 '.'A < sizes < 1 '.'A) are represented by thin points and 
large structures (1 '.'A < sizes < A f .'0) by thick ones. 

We have split the data sets into small (sizes < 1 '.'A) and large 
structures (1 '.'A < sizes < 4 "0). The results correspond to 
images filtered of p-modes. All time steps are included in these 
computations and their results are averaged. 

From these computations one can obtain the following infor- 
mation: 



• / con : For the family of large granular cells the curves show 
a sort of intensity level of saturation in both the centre of 
the granules and of the intergranular lanes. The rest of the 
transition can be described as a monotonic increase toward 
the granule. In case of small granular cells the general 
trend is essentially linear. The discontinuity (j um P) at the 
granular border (also visible for other line parameters) 
could be only apparent as a consequence of the uncertainty 
in the positions of the inflection points, which are not 
determined with sub-pixel precision. A glance at the upper 
panel of Fig. 12 reveals that the average of intensity values 
within each panel is not equal to 1 . This is caused by the 
fit of our spectra to the local continuum of the Liege- Atlas 
(Delbouille et al. 1973-1979) which do not correspond to 
the global continuum of this spectral region. 

• V: We find positive values (red shift) in the intergranular 
lanes and negative values (blue shift) in the granules as 
expected. The vertical speeds show increasing blue (red) 
shifts towards the centre of the granules (intergranular 
lanes). This is not in contradiction with findings of other 
authors indicating that upflows in large granules peak 
adjacent to the intergranular lanes; this refers to particular 
cases of very large granules whereas in our average of 
large granules sizes ranging from 1 '.'A to 4 "0 are included. 
Like in the case of 7 con the curves of Vqi for large granular 
cells present a saturation level of the vertical speed at the 
centre of the granules and the intergranular lanes and the 
strongest gradient is close to the granular borders. We 
observe an attenuation of the velocities with height in the 
photosphere that is in agreement with the literature and the 
results obtained before from the computation of response 
functions and correlation coefficients. 

• /: For the higher formed line (Line I) we find larger values 
of / in the intergranular lanes, confirming the reversal 
of temperature in higher layers. The line core intensities 
increase for both kinds of granular cells almost linearly 
from the centre of the granule to the central part of the 
intergranular lane. For Line III the difference of line core 
intensities between granules and intergranular lanes for 
small and large structures is very small, which confirms 
the previous results obtained by correlation and coherence 
analysis. We cannot expect large temperature fluctuations 
for line core intensities of this line, because most of the 
convective information related to Tni comes from layers 
around log T = —1.0, where the inversion of temperature 
occurs (see Table 2 and Fig. 8a). Note, however, that in 
Line III the global trend for the case of small granular cells 
is steeper than that for large granular ones, whereas in case 
of Line I the opposite occurs. 

• EW: The equivalent width (see Fig. 12) is larger in the 
granules as compared with the intergranular lanes. It 
shows a minimum located at the middle of the distance 
between the intergranular lane and the border, followed 
by a monotonic increase and a stabilisation towards the 
granular centre. Following Gadun et al. (1997) the steeper 
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temperature gradient in hot granular areas always produces 
deeper line profiles and respectively larger EW. 

• HW: A decrease of EW only caused by a change of tem- 
perature should also produce an decrease of H W. However, 
in Fig. 12 the opposite behaviour can be observed (small 
granular cells in Line I are an exception), a decrease of 
EW corresponds to an increase of HW. The simplest way 
to explain that is to include again the ingredients of turbu- 
lence or strong velocity gradients. Thus, the enhancement 
of HW can be interpreted as caused by an enhancement 
of turbulence from the border of the granule to the inter- 
granular lane. In case of Line I, the smoother increase of 
HW throughout the intergranular lane could be justified as 
a partial compensation by temperature changes, that is also 
reflected in a strong decrease of EW. Nevertheless, the ab- 
solute peak to peak variation (~ 8 mA) remains similar for 
both lines. These results are not in complete agreement with 
Nesis et al. (1993, 1999), who claim that in the deep photo- 
sphere, enhanced turbulence is concentrated predominantly 
near granular borders, while at higher layers the turbulence 
spreads out over the entire intergranular space. According 
to Solanki et al. (1996) the increase of HW can be justified 
as a consequence of the vertical velocity structure of the 
atmosphere in intergranular lanes, and turbulence must not 
necessarily be invoked. 

6. Detailed analysis of the variation of line 
parameters in space and time 

In this section detailed studies about the variation of line pa- 
rameters in space and time, separately for a small and a large 
granular cell (see Fig. 13) by averaging in sampled time or 
space, are represented. 
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Fig. 14. Convective component: Spatial scan across a small 
granular cell, where in the four lower panels Line III is repre- 
sented by solid lines, Line I by dashed ones. Granular borders 
are indicated by arrows. 
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Fig. 13. False colour image (left panel) and binary images 
of 

^con representing the selected areas for scans across small 
(small markers) and large granular cells (larger markers) in spa- 
tial (middle panel) and temporal (right panel) direction. 
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Fig. 15. Convective component: Spatial scan across a large 
granular cell, where in the four lower panels Line III is repre- 
sented by solid lines, Line I by dashed ones. Granular borders 
are indicated by arrows. 
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Figures 14 and 15 show the spatial variation of line 
parameters across a small (size < 1 "5) and a large (size ~ 
3 "0) granular cell, respectively. The borders of the granules 
are defined as the inflection points of I con , i.e. the extrema of 
the spatial gradient of intensities. We find the variation of line 
parameters more pronounced across the large granular cell. For 
V and HW this variation is attenuated with increasing height. 
Thus, these particular examples fit very well to the observed 
mean trend for V and HW in Fig. 12. 7 con and V show a broad 
maximum at the granular centre of the large cell. In lower 
layers the profile of the variation of vertical velocities in both 
cells shows a good (negative) correlation with the variation of 
/ con , whereas, as expected, this correspondence weakens for 
higher layers (see Table 4). In the larger structure, the minima 
of the variation of the line core intensity are shifted towards the 
granular borders. In concordance with our results in Sect. 5.2, 
higher (lower) values of equivalent width (full width at half 
maximum) above the granule (intergranular lane) are found. 

The temporal evolution of the several line parameters has 
also been studied for both cells (see Fig. 13). The maximum 
of the / con variation coincides in time with the extrema of the 
remaining line parameters. 

7. Conclusions 

The vertical structure of the solar photosphere has been 
investigated, analysing the global and spatial variation of line 
parameters obtained from a time series of ID spectra including 
two Fel lines with well separated formation heights. Thus 
we have been able to study how granules overshoot from the 
convection zone into the stable radiative atmosphere. For a 
qualitative interpretation of our results the formation heights 
of the several line parameters have been estimated, using the 
concept of the response function. The images, mapping the 
values of different line parameters in space and time have been 
filtered of p-modes. 

The convective motions of structures larger than 1 "0 pen- 
etrate up to the highest layers of the photosphere (~ 435 km). 
This conclusion stems from the significant correlation and co- 
herence between continuum intensity and vertical velocities for 
both lines. In case of the higher forming line (Line I) the corre- 
lation and the coherence decrease which can be interpreted as a 
less efficient penetration of convective motions into the high- 
est layers. In concordance, a significant correlation between 
the vertical velocities of both lines is found. This correlation 
is larger before filtering the p-modes, which implies a signifi- 
cant contribution of the oscillatory velocities at both layers. 

Regarding the intensity pattern, in very high layers (~ 435 
km) only structures larger than 2 "0 are connected with those at 
the continuum layer, but showing inverted brightness contrast. 
Thus the existence of an inversion of temperature is confirmed. 
We find such an inversion layer at ~ 140 km, coinciding with 
the height of formation of the line core intensity of Line III. 
This is derived from the fact that neither significant correlation 
nor coherence of the line core intensity (/ni) with continuum 
intensity (/ CO n) and line core velocity (Vqi) is found. 
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The phase shifts found between V\ and 7 con , /i, Vm could be 
interpreted as a possible evidence of gravity waves superposed 
on the convective motion (Deubner 1989) in higher layers. 

From the study of the dependence of line parameters to the 
distance from the granular borders, for both lines an enhance- 
ment of full width at half maximum (HW) towards the inter- 
granular lane is found. This goes together with a decrease of 
equivalent width (EW). A decrease of EW only caused by a 
change of temperature should produce a decrease of HW. The 
simplest way to explain our findings is to include the ingredi- 
ents of turbulence or strong velocity gradients in the intergran- 
ular lanes. 

For the continuum brightness and vertical velocities for 
the family of large granular cells (1 '.'A < sizes < 4 "0) a 
level of saturation in both, the centre of the granules and of 
the intergranular lanes is found. This reflects only an average 
behaviour and it does not necessarily correspond to cases of 
individual large granules, reported by several authors (e.g. 
Nesis et al. 1992, 1993, 1997; Hirzberger 2002) indicating that 
upflows in large granules peak adjacent to intergranular lanes. 

The study of the variation of line parameters across two 
particular granular cells, small and large, respectively, con- 
firms the mean behaviour, derived from the statistical analyses 
performed for the whole sample. 

Future aspects: In forthcoming papers, we will present 
studies of the spatial and temporal variation of physical 
parameters obtained from the inversion of our data using 
the inversion technique SIR of high resolution time series of 
spectroscopic data (Ruiz Cobo & del Toro Iniesta 1992, 1994). 
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